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C"| ' Abstract 
^ ' 

1 '— ^ ' This paper presents the perturbation theory for the double-sine-Gordon equation. We received the 

system of differential equations that shows the soliton parameters modification under perturbation's 
influence. In particular case A = the results of the research transform into well-known perturbation 
theory for the sine-Gordon equation. 
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Introduction 

^ ; The double-sine-Gordon (DSG) equation was studied in some physical problems. It 

describes spin waves in superfluid 3He, self-induced transparency in accounting degen- 
eracy of atomic levels P, electromagnetic waves propagation in semiconductor quantum 
CN ! supelattices 0, etc. The unperturbed DSG-equation can be written as 

O; ^-^ + sin^ + Asin2y. = 0. (1) 

\ DSG-equation can not be solved with the inverse scattering transform. But it is 

known that equation Q has the kink solution 

> : 

k*n ; y? = 2arctan ^V2A + 1 sinh" 1 (^=^=(x - x - ut)jj . (2) 

We shall call this solution as DSG-soliton, to differ it from the soliton of the sine- 
Gordon (SG) equation (particular case of the DSG-soliton, when A = 0). Interactions 
of solitary waves like (J2J) are not elastic [3J QJ 13] , they are accompanied by the radiation 
loss. The propagation media also influences the soliton parameters (its speed u and its 
location x ). This influence can be described as perturbation in (HJ). Thus, we have 

d 2 tp d 2 ip 

The solution (j2J) does not exactly satisfy (JSJ)- But if the perturbation is small we 
can take the DSG-soliton with slowly varying parameters as the solution. Then we find 
equations, which describe modifications of u and x . For this purpose we construct the 
perturbation theory. 
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Perturbation theory (main algorithm) 



This paper is not the first attempt to solve the perturbed DSG-equation. This problem 
was examined in assumption A < 1 [6 j using the inverse scattering transform for the 
perturbed sine-Gordon equation. But the paper did not show the DSG-soliton parameters 
modification. 

There is a number of methods for the perturbation theory constructing |7j- JTj. As it 
was shown in ^2] the perturbation theory for the SG-equation by McLaughlin and Scott 
describes the experimental results and the numerical solutions of the SG-equation better. 
Therefore, we use the algorithm represented in ^3]. For these purposes we modify and 
adopt this algorithm for the DSG-equation. 

Let's rewrite the DSG-equation in matrix form 
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We mark the column (ip, ip t ) as W 
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ft 



(4) 
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W is the function of the DSG-soliton parameters u and x . So, W = W(u, x ) or 
W = W(p), where p = (u,xq). We suppose that the solution of the equation (jHJ) can be 
presented as 

W = W + eW 1 . (6) 

Here Wo is the exact solution of (JIJ with the dependent on time parameters. 

If we insert this perturbed solution into (JIJ and linearize it, then we rewrite it in the 
following form 



L(Wx) = F(W ) 
Here L is the linear operator which is given as 



(7) 



L 



1 ° d+ ° _1 

1 * — d xx + cos (po + 2A cos 2</? 

The right part of (J7J) represents the 'effective force' 



F(W ) = /(Wo) 



1 2 dpi dW 



(8) 



(9) 



e fr{ dt dpi 

The boundedness condition for W\ is provided with the orthogonality F to the discrete 
subspace of the kernel Nd(L + ), where L + is the adjoint operator to L 



V 



1 
1 



d t + 



— d X x + cos tpo + 2A cos 2</? 
-1 



(10) 



Inserting the basis in Nd(L + ) and putting it into we derive the equation, that 
satisfies the orthogonality condition 



A f, dW \ dp, . 



:n) 



3=1 
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where (F(x), G(x)) = /f^ F T (x)G(x)dx. And b{ are the basis vectors of the adjoint space 



-1 

1 



dW 
dpi ' 



(12) 



The equation (|TT|) represents the perturbation theory for the DSG-soliton parameters. 

Analysis of the DSG-soliton motion 

And now let's consider the case of the DSG-soliton propagation. We have 

(13) 



(14) 
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ftu 



Thus, we receive the basis vectors of the adjoint space as 



—ftx 


, b 2 = 


-ftu 


fx 




fu 



It allows us rewriting the system (fTTj) 
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f(fo) 


f(fo) 



(15) 



The other summands in (fTTj) equal to zero. If we transfer these equations into integral 
form we receive 



/ (ftufx - ftxfu)dx) -jr = £ f f(ip )ip x dx, 

J — CO / OX J —CO 

dxo 



(16) 



(ftxfu - ftufx)dx 



dt 



f(<Po)tp u dx. 



These equations transform into a system of two equations. They describe the evolution 
of soliton parameters, its speed and location. 



du 
~dt ' 

dx 
~dt 



1 - n 2 2A+ 1 
~~ 2 k{\) 



fife 



cosh(6>) 



cosh'(#) + 2A 



u 



u — e- 



Vl-u 2 2A + l 
! jfe(A) 



f( V0 ) dx. 
J cosh 2 (6>) + 2A 



(x — Xq — ui), fc(A) = V2A + 1 H 1= tank 



-i 



2A 



Here 6 - ., ... ,. ,, , t v _. . . — ... , 

V 1 -u 2 V2X \ V 2A + 1 

If we put A = 0, then k(0) = 2 and this system reduces to a similar system for 

the perturbed SG-equation [T3]. It's easy to notice, that the parameters modification 

depends on the value of A and the type of perturbation function f(fo). This function 

can be written in different variations. But the most frequently appeared perturbation 

function contains three summands 
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a 



dtp 
dt 



/i[sin(y?) + \sxn(2<p)]S(x) 
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The influence of the constant energy pumping (the summand of 7), the energy losses 
(the summand with a) and the presence of the inhomogeneous region (like the microshort 
in the Josephson junctions) are represented in this perturbation function. If we take 
H = and substitute the expression of (jTSjl in (jT7jl . we shall get the following system of 
differential equations 



du 
~dt ' 

dx 
~~dT 



7T7 

2fc(A) 
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3/2 



au (l — w 



(19) 



u. 



Analyzing this system we derive the speed stability condition. To satisfy the condition 
the right part of the first equation must be equal to zero 



7T7 



(l-n 2 ) 3/2 -cm(l-w 2 ) 



2Jfe(A) 

So, we get the stabilized speed value 



0. 



(20) 
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1 + 



'2k(X)a 

. ^7 , 



(21) 



This expression shows that if the energy pumping is minimal (7 — > 0) the soliton 
speed reduces to 0. On the contrary if the DSG-soliton takes the energy much faster than 
it looses (7 ^> a) its speed increases up to 1. Putting A = we take the expression for the 
stabilized speed of the SG-soliton 13 . Any DSG-soliton reaches the speed of u under 
the action of energy pumping and energy loss without the difference of its initial speed. 

In case \i 7^ the system (|T7j) can be written as 



du 
~dt " 

dxQ 
~~d7 
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2k(X) 
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(22) 
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\iuxq v2A+T(7 ( 9o ) 
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(cosh 2 fx) +2A 



sinh(2x) 



Hereto W1 ^ 

If we put the values 7 = and a = 0, then we shall get the case of the pure 
interaction between the DSG-soliton and the inhomogeneity. The DSG-solitons with low 
speed reflects from the inhomogeneity. And high speed solitons overcomes this region. 
The results of interactions are shown on Figure ^ 

It is possible to find the condition of overcoming from the energy analysis. The 
Hamiltonian function of the unperturbed DSG-equation is defined as 



H 




cosy?) + -(1 



cos2y?) dx. 



(23) 
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Figure 1: The DSG-soliton interaction with the inhomogeneity. The parameters are A = 0.2 (left), 1.2 
(right); fj, = 0.5, 7 = a = 0. The initial soliton speed: 1) 0.2; 2) 0.4; 3) 0.45; 4) 0.5; 5) 0.6. 



Substituting (J2J) into (J2H)) we receive the expression for the full soliton energy 

4fc(A) 



E 



vT 



(24) 



To find the 'energy of rest' let's put the speed value which is equal to 0. Thus, we get 
the kinetic energy of the DSG-soliton 



E k = 4fc(A) 



(25) 



The inhomogeneous area can accumulate energy. This process can be counted as the 
correction to the DSG-Hamiltonian 



Hi = J \x \ {\ — cos (p) + — (1 — cos2(/?)J 8{x)dx. 
Having substituted (2) into the previous expression we get this correction 

' cosh #0 \ 



Hi = 2/j(2X + iy 



cosh 2 # + 2A, 



(26) 



(27) 



The maximal energy which can be accumulated by the inhomogeneity 

2fi, if - 0.5 < A < 0.5; 

(2A + 1) 2 



H„ 



(28) 



4A 



-fx, if A > 0.5. 



The value of the critical speed can be found from the condition of equality of the 
kinetic energy ()25|) and the energy accumulated in the inhomogeneity 
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Hmax + %H max k(\) 

\| [H max +Ak{\)Y 



(29) 



The one more thing must be mentioned about the DSG-soliton interaction with inho- 
mogeneity. If all media parameters (7, a, /x) are not equal to zero the phase path is not 
symmetric concerning the it-axis. This effect is shown on Figure 121 
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Figure 2: The DSG-soliton interaction with inhomogeneity in case a = 0.04, fi = 0.5, A = 1.2, 7: 1) 
0.030; 2) 0.032; 3) 0.035. 

This asymmetry is caused by the energy pumping presence. As soon as the soliton 
begins to loose its speed the equilibrium condition (|20|) breaks. From this moment the 
system has the positive balance of the energy gain. This process changes the DSG-soliton 
speed and, therefore, the phase path changes its shape. This effect takes place in the same 
problem for the SG-equation, but usually it isn't considered. The asymmetry is small if 
the media parameters a, 7 <C 1 or the soliton speed u — > 1. 

In addition to this effect we notice the soliton damping after reflection. The damping 
is caused by the influence of energy inflow that pulls the solitary wave to the positive 
direction along the x-axis. 

The numerical analysis shows that the view of the phase path depends on the value 
of the A-parameter (see Figure EJ) • 

It is necessary to notice that the condition ()28|) allows to classify the DSG-solitons. 
The DSG-soliton (exactly, its derivative) looks like one-humped solitary wave (like the 
SG-soliton), when A G (—0.5,0.5). If A > 0.5, the DSG-soliton is a two-humped wave. 
That's why the phase path has two turning points if . And the phase path has the only 
one extreme point in case A < 0.5. Thus, it is possible to evaluate the A-parameter 
observing the soliton behaviour. And vice versa. We can control the propagation process 
by changing the value of A. 
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Figure 3: The influence of A on the DSG-soliton motion in case 7 = 0.035, a — 0.04, fx — 0.5, A: 1) -0.2; 
2) 0.5; 3) 1.2; 4) 1.9. 

Conclusion 

In the conclusion we shall briefly observe some results. First of all, this paper contains 
the perturbation theory that has been adopted for the DSG-equation. The derived formu- 
las were used to consider the DSG-soliton propagation. It is shown that the DSG-soliton 
can be stabilized in the presence of the constant energy pumping and energy loss. We 
derive the condition of the stabilization. Also the DSG-soliton interaction with inhomo- 
geneity was studied. We notice two variants of interaction: reflection and overcoming. 
And so we receive the condition of their separation. The variant of the propagation 
depends on the parameters of the media: a, 7, /1 and A. 

In case A = all the results transform into the same expressions for the SG-equation 

[Si- 
Appendix A 

We represent the integrals that were used in the calculations. We use tanh -1 to mark 
the inverse hyperbolic tangent. To derive these expressions well-known integrals from [?] 
were taken. 
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(cosh 2 x + a) 2 
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The definite integrals were calculated. 

dx 1 



r 

j — i 



tanh 1 



cosh x + a ^ a ( a + 1) VV a + 1 



cosh x dx 1 1 / / a 

J-oo (cosh 2 x + a) 2 ~ a + 1 + y/E{a + l) 3 / 2 ^ \\ a + I 



<ix 1 , _i / / a \ 1 

■ tanh 



3 (cosh 2 x + a) 2 v/a(a + l) 3 / 2 \ y a + V a(a + 1) 

Appendix B 

There are the expressions for the function k(X): 



a) if A > 0, then fc(A) = V2A + 1 + -= tanh" 1 ' 



v/2A \V2A + 1/' 

b) if A = 0, then fc(A) = 2; 



1 / / — 2A 

c) if -0.5 < A < 0, then fc(A) = V2A + 1 + arctan 



V^2X IV2A + 1 
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